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Abstract. We write, for geometric index values, a probabilistic proof of the 
product formula for spherical Bessel functions. Our proof has the merit to 
carry over without any further effort to Bessel-type hypergeometric functions 
of one matrix argument. Moreover, the representative probability distribution 
involved in the matrix setting is shown to be closely related to matrix- variate 
normal distributions and to the symmetrization of upper-left corners of Haar- 
distributed orthogonal matrices. Once we did, we use the latter relation to 
perform a detailed analysis of this probability distribution. In case it is ab- 
solutely continuous with respect to Lebesgue measure on the space of real 
symmetric matrices, the product formula for Bessel-type hypergeometric func- 
tions of two matrix arguments is obtained from Weyl integration formula. 



1. Reminder and motivation 

The spherical Bessel function of index v is defined for all complex z and all 
> -1 by ([H]) 

^ (-1)' (z^^i 

where {v -f 1)/ := T{v -\- 1 + \)lT(y + 1) denotes the usual Pochhammer symbol. 
It provides a basic example of one-variable special function satisfying a product 
formula that opened the way to a rich harmonic analysis. More precisely, for 
V > —1/2 and nonnegative real numbers x,y, z, it is well known that 

(1-1) ji^{xy)j^{zy) = >(?yX,^(rfO> 

where r^^ is a compactly-supported probability distribution. Recall that for i/ > 
— 1/2, (jl.ip is a trivial consequence of the addition Theorem for Bessel functions 
(see for instance Chapter XI in [16]) while it obviously holds for = —1/2 since 
j-i/2{z) = cos(z). Nevertheless, for integer p > 1 and for the so-called geometrical 
index values v — {p/2) — 1, (jl.ip may be derived from the following Poisson-type 
integral representation 

(1.2) J(p/2)-i(IH) = / e^^^'^^a.ids), veRP, 
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where cri is the uniform distribution on the unit sphere S^^^ and (•,•), |-| are 
respectively the Euchdean inner product and the associated Euchdean norm in Rp. 
Indeed, if we set \v\ = y, then 



where cFx^az are the uniform distributions on spheres of radii x,z respectively. 
But according to pjj Corollary 5.2 p. 1149, the probability distribution (t^ -ka^ is 
absolutely continuous with respect to the Lebesgue measure in W and due to its 
rotational invariance it has a radial density. The use of spherical coordinates yields 
then (jl.ip . Avoiding techniques from differential geometry like the ones used to 
prove the absolute continuity of we write a probabilistic proof of (II. ip for 

geometric index values and supply a probabilistic interpretation of ri^/^^* ^ . Our 
starting point is the elementary fact that the conditional distribution of a standard 
normal vector N in W given its radius |A^| is the uniform distribution on the 
sphere of radius |iV|. The product of two spherical Bessel functions turns towards 
the conditional independence of two independent standard normal vectors N\ , Ni 
relative to the ct- field generated by their radii |A^i|, |7V2| ([13]). The representative 
probability distribution ri^/^^"^ is then seen to be the conditional distribution of 
the radial part |7Vi + A^2| given (|iVi| = x, |7V2| = z). In fact, A^i + A^2 is again 
distributed as a standard Gaussian vector (up to a constant) and its angular part is 
independent from both radii |A^i| and |A^2|- The reader will easily realize from the 
ingredients needed in the proof that choosing any multivariate stable distribution 
in WP whose density is a radial function does not alter our proof. But the Fourier 
transform of a radial function is again radial therefore the choice restricts uniquely 
to isotropic or rotationally invariant stable distributions (whose Levy exponents are 
given up to a constant by w |?;|", a G (0, 2], [IS] p. 86). 

Our proof has also the merit to carry over after mild modifications to some 
matrix analogues of spherical Bessel functions. Those we consider here are known 
as Bessel-type hypergeometric functions of one and two m x m real symmetric 
matrix arguments. The product formulas we obtain are valid for geometrical index 
values and are those derived in [T^ using hypergroup theory, in the particular case 
of the real division algebra. This is by no means a loss of generality since product 
formulas over the division algebra C may be easily derived along the same lines. For 
functions of one matrix argument, the proof is identical to that written for j(p/2)-i. 
Besides, the representative probability distribution is seen to be the conditional 
distribution of the radial part of the sum of two independent pxm{p> m) standard 
matrix-variate normal distributions given the radial part of each. We shall prove 
that this conditional distribution is closely related to the distribution of the mx m 
upper-left corner of an orthogonal matrix of size p, whence its absolute continuity 
(with respect to Lebesgue measure) is deduced for p > m + 1. For these values 
of p, one easily derives the product formula for functions of two arguments using 
Weyl integration formula for the space of real symmetric matrices. As a matter 
of fact, the corresponding representative probability distribution has an analogous 
description in terms of singular values rather than matrices. Besides, when p > 2m, 
a result due to B. Collins provides a detailed description of the distribution of the 
upper-left corner of an orthogonal matrix, agreeing with the variable change formula 
given in Lemma 3.7 p. 495 in [9] and reproved in Corollary 3.3 p. 762 in [M]. Note 
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finally that since Bessel-type hypergeometric functions of two matrix arguments we 
consider here are instances of generalized Bessel functions associated with _B-type 
root systems, then our approach resembles the one carried for proving Theorem 
5.16 (ii) in il . 

The paper is organized as follows. In the next section, we consider spherical 
Bessel functions j(p/2)-i and prove for geometric index values. In section 3, 
we extend our proof to Bessel-type hypergeometric functions of one real symmet- 
ric matrix argument. In the last section, we perform a detailed analysis of the 
representative probability distribution: it is absolutely continuous for p > m + 1 
and its density enjoys a certain averaged bi-invariance property with respect to the 
orthogonal group. The product formula for functions of two real symmetric matrix 
arguments follows from Weyl integration formula. 



2. Product formula for spherical Bessel functions 

All random variables occuring below are defined on some probability space 
{il,J^,F) and we denote E the corresponding expectation. Furthermore, for the 
cr-field cr{X) generated by a random variable X, we write 

E[-\X] for E[-\a{X)], 

and we recall that all equalities involving conditional expectations hold P-almost 
surely. Let iV be a standard normal vectoiQ in W and let N = RQ be its polar 
decomposition (i? > and 8 G S^^^). Then, R and Q are independent and 8 is 
uniformly distributed on S^^^. It follows that for any t; e 



E 



J{v,Rs) 



CTl(ds) = j(p/2)-li\v\R)- 



In fact, if X, Y are independent random variables valued in some measurable spaces 
and if T>y stands for the distribution of Y, then 

E[f{X,Y)\X] = I f{X,y)VY{dy) 

for any bounded Borel function / (see J,3J p. 108 Exercice 4.27). 

Now, let Ni , N2 be two independent standard normal vectors in Rp with polar 
decompositions Ni = i?i8i,iV2 = i?282 respectively, and consider the product a- 
field a{Ri, R2) generated by i?2- Then, the independence of A'^i and N2 implies 
that ([13J) 



E 




Ri 


= E 




i?2 


E 




R2 


^ E 




i?2 



Besides, Ni,N2 are conditionally independent relative to o-(i?i,i?2) (see [T3j p. 109 
Exercice 4.32). In fact, one has for any bounded Borel function / : Rp — > M 



Thus 



IE[/(A^2)|iVi,i?i,i?2 



E 



f{N2)\R2 

g^(i,,iV2)|^ J =E 



E 



/(7V2)|i?i,i?2 



Its coordinates are independent centered normal distributions with unit variance. 
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Write A^i + N2 RsQa, then A^i + N2 is (up to a constant factor) a standard 
normal vector so that Q3 is uniformly distributed on S^^^ and is independent from 
R3. We claim that: 

Proposition 2.1. O3 is independent from cr(i?i,i?2)- 



Proof. Let / : S-P-i ^ R,5 : K 
independence of A^i , N2 yields 



be bounded Borel functions, then the 



,[/(e3)5(i?i,i?2)] =E 



/ 



NI + N2 



00 /'OO 



Jo 



\Ni+N2\ 
F{ri,r2)dridr2 



gilNiim) 



/ 



5p-1xSp-i 



|ri6'i +r26l2| 



ai(d0i)ai(rf02), 



where 



F(ri,r2) := {rir2Y-^e-^-"+-"^/^ g{r,,r2). 

Let 

t^n ,r2 (^^) be the pushforward of (T\ ® (j\ under the map 

riOi + r202 



|ri0i+r202r 



then 



L 



f 



al(^^0l)al(^^^^2) 



5P-1 



But h'ri,r2 is obviously invariant under the action of 0(p), therefore i'ri,r2 = ""i since 
(Ti is the unique distribution on Sp~^ enjoying the rotational invariance property. 

□ 

We also need the following lemma: 

Lemma 2.2. Let V,X,Y be random variables such that Y and {X,V) are inde- 
pendent. Then, for any bounded Borel function f 

E[f{X,Y)\V] = I E[f{X,y)\V]Vr{dy). 

Proof. This fact is easily proved for bounded functions f(x,y) = g{x)h(]j) and 
then extended to bounded Borel functions using the monotone class Theorem ([T^ 
p.5). □ 

Combining the proposition and the lemma, one gets 



E 



E 



5P-1 



Finally, let M-RsK-Ri.iig) be a regular version of the conditional distribution of R^ 
given {Ri,R2), then Fubini Theorem entails 



Thus, (jl.ll) is proved and T; 



(p/2)-l 



as explained in the following remark. 



fits t''Ri\{B.i,R2) tii6 event = x,R2 = z} 
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Remark 2.3. Let <E> be the angle between 61, 82: cos$ = (81, 82). Then 

i?3 = 



Ri 



RI + 2R1R2 cos$. 



But the independence of 81, 82 entails for any real w 



E[e 



iw cos 



e™<"'*Vi(ds)CTi(di) 

gp-l J5P-I 

jip/2)-i{w\t\)<^i{d.t) 

r(p/2) 



SP 



i(p/2)-l(w) 



r(i/2)r((p-i)/2) 

where we used Lemma 5.4.4 p. 195 in 5 . Performing the variable change 

one recovers the density of ri^/^-* ^ derived in Proposition A. 5. p. 1153 in |llj . 

3. Product formula for Bessel-type hypergeometric functions of 
one real symmetric matrix argument 

In this section, we consider matrix- variate normal distributions rather than vec- 
tors. Doing so leads to a product formula for Bessel-type hypergeometric fimctions 
of one real symmetric matrix argument (see below). To this end, we recall from ^ 
Ch.L the following needed facts. Let p > ni > 1 and let be a real matrix-variate 
p X m standard normal distribution, that is a p x m matrix whose entries are inde- 
pendent centered normal distributions with unit variance. Then N admits almost 
surely a unique polar decomposition N = Z{N'^ NY^"^ := ZH. Moreover, Z and 
H are independent, H is almost surely invertible and Z is uniformly distributed on 
the real Stiefel manifold 



-'p^m 



where M, 



is the space of p x to real matrices. Let 0{p) be the orthogonal 
group, then Ep_m is a homogeneous space Ep „i « 0(p) / 0{p-'m). It thereby admits 
a unique 0(p)-invariant distribution we shall denote Up^m- More precisely, Up^m is 
the pushforward of the Haar distribution on 0{p) under the map 



O ^ Oep 
Hence, for any C G Afp.„j(M) 



'0„ 



ids) 



ids). 



Now, let iVi , A'^2 be two independent p x m matrix-variate standard normal distri- 
butions with corresponding polar decomposition Ni = ZiHi,N2 = 2^2-^2- Then, 
by considering the product (T-field cr(_ff i,_ff2) generated by Hi,H2 we easily derive 



(3.1) 



E 



Since Ni + N2 is up to a constant factor a p x m matrix-variate standard normal 
distribution, then it admits almost surely a polar decomposition A^i + iV2 = ^siJs, 
where Z3 is uniformly distributed on Ep^m and is independent from H3. Similarly 
to the case m — 1, one proves that Z3 is independent from a{Hi,H2) (analogue of 
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proposition I2.ip using the foUowing variable change formula ([7], Prop. XVI. 2.1. 
p. 351): let dA be the Lebesgue measure on A/p^m(R), let S'+(R) be the set of real 
positive definite matrices with Lebesgue measure dr and 7 = (15/2) — 1 — [m{m — 
l)]/2. Then 



f{A)dA 



Mp,„(R) 



f{sy/r) [det(r)]''(Tp,„i {ds)dr. 



Accordingly and with the help of lemma 12. 2[ one gets 



and if P^H3\(Hi,h.2) is the conditional distribution of H-^ given {Hi, H2), then Fubini 
Theorem entails 



E 



Using [9|, (3.5) p. 493, one sees that 



E 



where qFi is the Bessel-type hypergeometric function of one real symmetric argu- 
ment and of geometrical index value (p/2) (it reduces when to = 1 to i(p/2)-i, ,10.). 
Finally, p.ip yields the product formula 



oFi —HiC^ CHij qFi —H2C^CH- 



Now, we proceed to 



4. Absolute continuity of M-f/3|(-f/i,ff2) and Product formula for 
Bessel-type hypergeometric functions of two matrix arguments 

4.1. Absolute continuity of /iff, 1(^1,^2)- contrast to the case to = 1, the 
absolute-continuity of ^J^H3\{Hl,H2) not obvious and needs a careful analysis we 
perform below: 

Proposition 4.1. For any p > to -I- 1, fJ-H3\{H2,Hi) absolutely continuous with 
respect to the Lebsegue measure on 5'm(M) and its density, say f(^jj^ [{^•■^{A), satisfies: 
(4.1) 

f{o^H,oT.02H20^)iOlA03)dO<S)dO= I 

0(m)xO(m) JO{m)xO(m) 

almost surely for any O3 G 0{m), where dO is the Haar distribution on 0{m). For 
p = m, it is singular. 



f{O^H^Oj ,02H20'r){A)dO®dO 



Proof. Since 

(ila)' = {Hif + {H2f + HiZlZ2H2 + H2Z^ ZiHi 
then t^H3\(H2,Hi) is the pushforward of Cp,™ ® crp,m under the map 



(Zi, Z2) ^ V (^1)' + (^2)2 + HiZ( Z2H2 + H2Zi ZiHi 
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for fixed Hi,H2, where for a positive semi-definite matrix A, \fA is its square root. 
But from the very definition of CTp.m, y^H^\{Hx,H:i) is the pushforward of the Haar 
distribution dO ® dO on 0{p) x 0{p) under the map 

or equivalently 

since dO is invariant under O i~> O-'". Besides, the random variable O1O2 G is 
Haar distributed since it is 0(p) -invariant. As a matter of fact, A* 1(^1,^2) the 
pushforward of dO under the map 

O ^ ^{H^Y + {H2f + H^el^0ep^mH2 + H2el^0^ep,rnHi. 
Now observe that for fixed Hi,H2, 

O ^ {Hif + {H2f + Hiel,^Oep^„M2 + H2el^O'^ Cp^^Mi 

is a affine map from 0{p) into 5m(R), therefore is hpschitzian whose differential 
map is constant. Moreover 0{p) and 5'm(K) are real analytic manifolds such that 
dim 0{p) — p{p — l)/2, dim S'm(M) — m{m + l)/2. As a matter of fact 

• If p = m + 1, then dim 0{m + 1) = chm S'm(M) and Theorem 3.2.5 p. 244 
in [5] implies that the pushforward of the Haar distribution on 0{p) under 
this map is absolutely continuous with respect to the Lebesgue measure on 
5™(R). 

• lip>m + 2, then dim 0{p) > dim 5„i(R) and Theorem 3.2.12 p. 249 in 
[5] yields the same conclusion. 

Now, for any d, 02,03 € 0(m), f(OiHiOj ,02H20^}{03 AO3) is the density of the 
random variable (for fixed Hi, H2) 

OMHifO^Ol + OMH2?OlOl+ 

0:iOiHiOl Zj Z2O2H2OI0I + 0302H20jZ^ZiOiHiO'[Oj 
which can be written as 

(030i)(i/i)2(Of Oj) + (0302)(i/2)'(02^03^) + 

(030i)i/i(OfOj)(Zi03^)^(Z203^)(0302)iJ2(O^Oj) 

+ (0302)1/2(02^ Oj)(Z20[)^(ZiO^)(030i)iIi(Of 03^). 

But since ap^m is invariant under the right action of 0{m) ([2j p. 28) and since the 
Haar distribution dO is 0(m)-bi- invariant, then the f[Hi,H2) satisfies (|4.1I) . Finally, 
since dim 0{m) < dim S'm(R) then Theorem 3.2.5 in [51 shows that for p — m, 
P^H-i\(Hi,H2) is singular with respect to the Lebesgue measure on S',„(]R). □ 

Remark 4.2. Note that 

^p.ra^^P,^ — A^ Op— m,p— m 

where A^ is the upper-left m x m corner of the orthogonal matrix O. According 
to [1], Remark 2.1. p. 118, if p > 2m then the distribution of Am is absolutely 
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continuous with respect to the Lebesgue measure on Mm m{R)- its density is given 

by 

det(I™-AA^)(f-2"-i)/2l{ll^ll<ij 

where 1 1 • 1 1 is the matrix norm induced by the Euclidian norm | • | . This fact should 
be compared with Lemma 3.7 p. 495 in [5]. 

4.2. Product formula for functions of two matrix arguments. Let p > m + 1 

so that I^H3\{H2,Hi) is absolutely continuous with respect to Lebesgue measure on 
5'm(R)- Then one derives a product formula for the Bessel-type hypergeometric 
functions of two real symmetric matrix arguments and of geometrical index values 
p/2,p > 1: if A is a real positive semi-definite matrix and C S Mp^„i(R), then these 
functions are related to those of one real symmetric matrix argument by 

(4.2) oFi (1;A--C^C) ^ [ oF, -OVAO^ {C^C)0VA0^) dO 

' Jo{m) ' 

where dO is now the Haar distribution on 0(m) (Theorem 7.3.3 p. 260 in 10 ). 
Keeping the same notations used in the previous section, one has 



^2 / Jo(rn) L 



dO 



loir 

which in turn implies that for any positive semi-definite matrices A, B and any 
CGAfp,„(R) 



oi^i (|; A; -C^C) oFi (|; i?; ~C^c) = 

JO(m)xO{m) JS+{R) / i' 2' 

Recall now that f(Hi,H2) denote the density of P'H3\{Hi,H2)- Then Weyl integration 
formula for ^^(R) ([Bl Theorem 10.1.1. p.232), and Fubini Theorem entail 



I (f ' -^C^^^) fi.O,^Oj.02VBO-)im®dO®dO = Cm f 

J0(m)x0{m) JS+{R) ^ Jo{m)xO{m) 

Fi (^;~ODO'^iC^C)ODO^) f(o,VAo^ 02^o^)^ODO^)V{D)dD®dOmO®dO 
= cm I oFi -ODO'^{C^C)ODO'' 



/ fio.^oT,02^o-)iD) dO®do\ V{D)dD ® dO 

JoMxOhn) 1 ' ^ 2 I 



lo(m)xO{m) 1 ' ^ 2 I J 

where D = diag(Ai > A2 > . . . > A„i) is a positive definite diagonal matrix 

m 

V{D):= II {Xn-Xj), dD^Y[dX,, 

1 < n < J < m j — ^ 

and Cm is a normalizing constant. By the virtue of (14. 2p . one gets 



oFi (|; A; -C^C) oFi (|; i?; ^C^C 



, f oi^i (^;D^;-C^C) KA.,B{D)dD 
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where 

KA,b{D) := V^(£»)l|Ai>...>A„>0} / f(o,VAor,02VBO?)i^)^0 ^ dO. 

Joim)xO{m) ^ ^ ^ 

Finally, one performs a change of variable Aj ^ \/Ai, 1 < i < m in order to get the 
product formula: 

oFi(|;^;-C^c) ,F,[^;B;-C^c) = 

2'"A,>...>A„>o V2' ) VAi...A,„fJ 
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